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Abstract 



We show that renormaUzation group improvement, making use of analyt- 
icity and the structure of the operator product expansion, combined with a 
non-perturbative expansion method allows one to describe quarkonium states 
via QCD sum rules without the need for explicit power corrections. This 
technique also allows an accurate determination of parameters related to the 
inclusive decay of the r lepton. 
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While the QCD sum rules technique |T|] provides a powerful method for extracting 
hadronic characteristics from QCD, there is the drawback that a description of the large 
distance behaviour via power corrections associated with quark and gluon condensates in- 
troduces extra phenomenological parameters. In this Letter we shall argue that the non- 
perturbative expansion technique of Ref. 0, based on a new small expansion parameter, 
combined with a particular renormalization group improvement which respects both ana- 
lyticity and the structure of the operator product expansion, enables us to describe some 
fundamental hadronic properties without the explicit introduction of condensates. As a 
consequence of the retention of the analytic properties of the running coupling, which is 
essentially equivalent to introducing a renormalon contribution, one may also evaluate the 
Rr ratio for the inclusive semi-leptonic decay of the r-lepton both by integration over the 
physical region and via use of Cauchy's theorem 0. 

The non-perturbative technique we use involves the so-called "floating" or "variational" 
series. In quantum mechanics this series has been applied successfully to the anharmonic 
oscillator @,|^. The variational series also appears in the linear (5-expansion |P and in 
variational perturbation theory ||^,|^ . In the case of zero- and one- dimensional field theories 
there exist rigorous proofs of the convergence of such series [§,0. When one considers 
QCD, the same ideas lead to an expansion for the generating functional in terms of a new 
parameter a, related to the coupling by p[ 



9^ ^ 1 a" 
{A-nf C{l-a] 



A-7f^ = ^77^, (1) 



where C is a positive constant which can be found from meson spectroscopy. It is clear that 
for all values of the coupling A > the expansion parameter a lies in the range < a < 1. 

In order to illustrate our technique, let us consider the Adler D-function D{Q^) = 
—Q'^dU/dQ^ corresponding to the vector hadronic correlator in the massless case. The 
two-loop perturbative approximation is given by D{t,X) = 1 + 4A(/i2), where t = QV/^^- 
Standard renormalization group (RG) improvement leads to the substitution A(/i^) — *• 
A(t, A), which implies a summation of the leading logarithmic contributions. However, due 
to the Landau pole of the running coupling at = ^^qcd this substitution breaks the 
analytic properties of the Z)-function in the complex = —Q^ plane, namely that the 
D-function should only have a cut on the positive real axis. 

We may correct this feature by noting that the above solution of the renormalization 
group equation is not unique. The general solution is a function of the running coupling 
with the asymptotic behaviour 1 + 4A, for small A. To maintain the analytic properties of 
the D-function we can write it as the dispersion integral of R{s) = (l/7r)Imn(s + ie), namely 

and use RG improvement on the integrand rather than D itself. This method leads to 
D{t, A) = 1 + 4Aeff(t, A), where, with r = s/Q"^, 

dr A(t, A) 

Jo (l + r)2l + A(t,A)/3olnr' ^ ^ 

This has the Borel representation 
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/■oo — 

Aeff(t,A) = / d6e-WA)5(5) ^ (4) 
Jo 

with 

B{b) = r(l + 6/5o)r(l - b(3o) . (5) 

Here f3o = 11 — 2/3 A^/ is the first coefficient of the /3-function, and A'^y is the number of 
active flavours. Thus, in the Borel plane there are singularities at 6/3o = —1,-2,... and 
6/?o = 1)2, ... corresponding to ultraviolet and infrared (IR) renormalons respectively. 

The ffist IR singularity at 6/5o = 1 is probably absent since there is no corresponding 
operator in the operator product expansion. Although this issue is not currently settled, it 
seems reasonable to assume that the ffist IR renormalon occurs at 6 = 2//5o, and we would 
like to use this property of the operator product expansion as an additional constraint on the 
choice of solution to the renormalization group equation. This can be achieved by integrating 
by parts in Eq. (]^), using the fact that to two-loop order R{s) is a constant, and applying 
the same RG improvement to the new integrand, to obtain the following expression for Aefr: 

Aefr(t,A) = r druir) J'l^'^^l , (6) 
Jo 1 + X{kt, Ajpo Inr 

in which the factor k reflects the renormalization scheme ambiguity and the function 

coir) = (7) 

describes the distribution of virtuality usually associated with renormalon chains. The 
function coincides with the function used in [|1^ and is numerically very close to that 



found in [|l^]. The Borel transform of (P) has the form 

B{b) = r(l + bf3o)T{2 - b[3o) . (8) 

Thus in this representation for X^s the positions of all ultraviolet singularities remain un- 
changed, but the ffist IR renormalon singularity at 6 = l//3o is absent. 

In summary, a representation for the effective coupling, and consequently the D-function, 
which manifests renormalon-type characteristics can be obtained as a particular RG im- 
provement of the lowest order radiative corrections which takes into account both the an- 
alytic properties and the structure of the operator product expansion. In order to render 
Eq. d^) integrable we must combine this method with the non-perturbative a-expansion of 
Ref. [0] in which from the beginning the running coupling has no ghost pole. Effectively, 
the representation for the D-function obtained in such a way coincides with a technique 



explicitly introducing power corrections |T3, and we can, in principle, describe hadronic 



parameters using, say, the method of QCD sum rules. 

However, retention of the correct analytic properties under RG improvement also allows 
for the possibility of considering inclusive r decay, and we shall consider this ffist. Using 
RG improvement following the procedure described above in the context of the a-expansion 
leads to 

K^{q') = -q' da \,, X{a), (9) 
Jo [cr — q'^Y 
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where to leading order 



A 



C 



[\ + 3a) 



(10) 



The renormahzation group defines the running of a(cT) as the solution of the transcendental 
equation P] 



a = Gq exp 



C 

2Po 



if (a) - /(ao) 



(11) 



where 



^ 2 6 18 1 

/(a) = 481na- — 

a 11 1 — a 

624 , , 5184 , 9 , 

H In (1 - a) H n (1 + - a) 

121 ^ ^121 ^ 2 ^ 



(12) 



The parameter ao and the virtuality ctq in Eq. ([IT]) are defined by some renormahzation 
point for the effective coupling. 

We now apply the representation (||) to inclusive r-decay. The starting point is the 
expression 



Rr = 2 

where R{s) = R^R{s) and 



ds 



2s 



77^ 1-77^ l + 77^^(^) 



M2 



M2 



M2 



(13) 



/?° = 3(|Kd|' + |Ks|')^w (14) 

in which the electroweak factor Sew = 1-0194 and the CKM matrix elements are \Vud\ = 
0.9753, \Vus\ = 0.221, taken from Ref. |]15[. Then, in order to isolate the QCD correction to 
Rr, we write Rr = i?° (1 + ARr). 

The effective coupling is an analytic function of s in the complex s-plane with a cut 
along the positive real axis. AR^. may be written as the contour integral 



A ^ di f dz ,^ 
AR^ = — <p — (1 
27ri J\z\=i z 



(15) 



where cii = 4 is the 2-loop coefficient of the D-function. 

Substituting Eq. (|) into Eq. ([TSl) and using Cauchy's theorem we obtain 



AR, = 12d, ^{^\ il- ^\ Mki 



(16) 



in which the factor k again parametrizes the renormahzation scheme. In what follows we 
shall always use the MS scheme, in which [|T^ k = exp(— 5/3). Note that in comparison 
with ([T3|), use of the representation (|^) for the coupling modifies the kinematic factor so 
that the maximum now occurs near s = (2/3)M^. 
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Taking as input the experimental value of R^^"^ = 3.56 ± 0.03 |T^, three active quark 
flavours and the variational parameter C = 4.1 as in we find 

a,(M^) = 4:71 ~X{M^) = 0.339 ± 0.015. (17) 

which differs significantly from that obtained (^^(M^) = 0.40 in leading order 0) without 
the renormalon-inspired representation (|^) for the coupling. The new method, applying the 
matching procedure in the physical s-channel and using standard heavy quark masses pO 
leads to Rz = 20.90 ± 0.03, which agrees well with experimental data 



As mentioned above, the requirement of analyticity naturally generates power correc- 
tions. Consequently, it is natural to investigate whether this technique allows us to describe 
meson parameters using the QCD sum rules approach. Consider, for example, the hadronic 
correlator n(s) corresponding to the vector current [|l^]. The imaginary part of n(s) to 



order in the non-perturbative expansion of Ref. is given by 

1 



Imn(s) = — 
47r 



nW(s) + 4An«(s) , (18) 



where A is given by Eq. (|T0|) and for fl^") and II*-^-* we use the perturbative two-loop expres- 
sions given in [|19|. It should be stressed that while the a-expansion technique makes use of 



the perturbative coefficients, the structure of the new expansion is fundamentally different. 
That it is indeed valid to use perturbative formulae in a non-perturbative calculation can 
be seen by noting that for the cc bound states to be considered, the dominant contribution 
to the moment integrals comes from an energy scale well above Aqcd- 
The initial expression for the first power moment is given by 



dQ 

I rco Imn(M 

" ^ io ''(Q2 + a + 4m2) 



da ,^, , , (19) 



where u = ^{(t/{(t + Am?)). In analogy with the massless case we shall associate the 
parameter cr = s — Am? with virtuality. The zeroth-order contribution is given by 



CO 



Mf\Q') = -^J da ^^^ ^^'^^'] (20) 
^ ' Att^ Jo (Q^ + o- + 4m2)2' ^ ^ 

while the order 0(A) correction term associated with 11*^^^ can be rewritten, again in analogy 
with the massless 

Mi"iQ^) = \r'''' J"^*'f\.M u), (21) 



Jo (Q^ + o- + 4m2) 

where 



/•u Oiij 

^{u) = il-u')l dv--^^U^'\v). (22) 
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The structure of the integrand in the expression for Mi allows us to simulate a summation of 
all terms by placing the term proportional to A in the denominator (cf. the mass correction 
for the propagator). As a result we have the expression 

1 /-oo T\(^)(v] 

-i^j/o '^w^wkr-' '''' 

where, after RG improvement, W{a) has the form 



W{a) = {a + 4m\k(T)) 



AX{ka) 



n(o)i 



u 



(24) 



and now u = {a / {a + Am'^ {ka))) . Note that although this procedure introduces correction 
terms of order A^, it is valid a posteriori, as the dominant contribution to the moment 
integral comes from the region where A is small |^ . 

Evaluation of this expression uses the running expansion parameter a(cr) obtained from 
Eq. (jn\), while the running mass, fn{a) = ■mQ(X(a) / X(ao))"'°^^° , has the standard perturba- 
tive form, with 70 the first coefficient of the anomalous dimension. However, in this case 
A (a) is expressed via Eq. (|l|) in terms of the running expansion parameter a. The moments 
of the vector correlator for general n are then given by 

A crucial feature of this approach is that the function W{a) achieves a minimum at some 
a = a, and for large n the dominant contribution comes from this point. Defining the 
ratio of the moments, RniQ"^) = M„_i(Q^)/M„((5^), we see that for large n this ratio tends 
asymptotically to its saddle-point approximation, i.e. 

K(g')™g' + w^(^). (26) 

On the other hand the corresponding phenomenological ratio for mesons, R^^'^{Q'^), has the 
form + at large n, where M is the mass of the first resonance in the relevant channel. 
Consequently, we obtain 



M = ^JW{a). (27) 

In [jl]] the moments were considered at = 0, while the case of 7^ has also been 
considered in the literature (see [jl9|| for a review). It is an important feature of this approach 



that the analysis is independent of Q^, as the functional form of the Q^-dependence for large 
n is the same for both the QCD and phenomenological moment ratios. 

Estimates for the cc bound state masses for the vector and axial vector channels, corre- 
sponding to (^) are shown in Fig. |l] for values of a near a. The analysis has been presented 
thus far for the vector channel only. The calculation in the axial-vector case is similar, the 
modifications to the moment expressions being standard, and will not be discussed here (see 
e.g. fT^)- Keeping C = 4.1 and taking ao at the r mass scale from Eq. (|T^ we have only 
one free parameter, the quark mass mo. With mo = 1.483 GeV we obtain a good fit to the 
experimental masses of both J/ip{lS) and Xci(l-P), as shown in Fig. |l|. 
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FIG. 1. A plot of ^/W{a) versus a for the vector (solid curve), and axial- vector (dashed curve) 
currents, the minimum being the asymptotic limit of the moment ratios for large n. For comparison, 
the straight lines are the corresponding experimental cc bound state masses [pO[. 



Note that the parameter C = 4.1 used here was found from a fit to the /3-function and 
corresponds to the linear growth of the quark-antiquark potential. The present calculation 
is consistent with this derivation, since the dominant contribution to the moment integrals 
arises from a region where the system is still reasonably non-relativistic (m^ ~ 0.3), 

In conclusion, we have presented in this Letter a technique for obtaining quarkonium and 
r-decay parameters which makes use of RG improvement and a non-perturbative expansion 
which removes the Landau pole in the running coupling. It appears that careful control of RG 
improvement, to ensure the correct analytic properties, naturally induces power corrections 
|l^ required for the description of meson parameters without the need for their explicit 



introduction. In the case of QCD sum rules the fact that this technique involves only three 
parameters, two being fundamental to QCD, i.e. the coupling constant ao and the quark 
mass mo defined at some energy scale, suggests that it is applicable to the study of other 
channels in the cc family and also to consideration of other heavy quarkonium states. Such 



an extension is currently under investigation and the results will be published elsewhere p2 
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